^/-statistics of Ornstein-Uhlenbeck branching particle system 
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Abstract 



We consider a branching particle system consisting of particles moving according to the Ornstein- 
Uhlenbeck process in W'- and undergoing a binary, supercritical branching with a constant rate A > 0. 
p I , This system is known to fulfil a law of large numbers (under exponential scaling) . Recently the question 

, of the corresponding central limit theorem has been addressed. It turns out that the normalization and 

, form of the limit in the CLT fall into three qualitatively different regimes, depending on the relation 

between the branching intensity and the parameters of the Orstein-Uhlenbeck process. In the present 
paper we extend those results to U -statistics of the system proving a law of large numbers and a central 
' limit theorem. 
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1 Introduction 



We consider a single particle located at time t = at x £ M , moving according to the Orstein-Uhlenbeck 
process and branching after exponential time independent of the spatial movement. The branching is 
supercritical and given by the generating function 

F{s) :=ps'^ + {l-p), P>^- 

The offspring particles follow the same dynamics (independently of each other). We will refer to this 
system of particles as the OU branching process and denote it by X = {X(}^>q. 

Formally, we identify the system with the empirical process, i.e. X takes values in the space of Borel 
measures on M*^ and for each Borel set A, Xt{A) is the (random) number of particles at time t in A. 
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for Polish Science. 
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It is well known (see e.g. |15| ) that the system satisfies the law of large numbers, i.e. for all bounded 
continuous functions, conditionally on the set of non-extinction 



Xtruxt,f)^{^,f) 



a.s. 



(1) 



where \Xt\ is the number of particles at time t, {Xt{l), Xt{2), . . . ,Xt{\Xt\)} are their positions, {Xt, f) := 
X^i=i fi-^tii)) and is the invariant measure of the Orstein-Uhlenbeck process. 

In a recent article [T], we investigated second order behaviour of this system and proved central limit 
theorems corresponding to ([1]). It turns out that the behaviour of the system may fall into three qualita- 
tively different categories, depending on the relation between the branching intensity and the parameters 
of the Orstein-Uhlenbeck process. 

In the present article we extend these results on the LLN and CLT to the case of [/-statistics of the 
system of arbitrary order n > 1, i.e. to random variables of the form 



(note that with this notation {Xf, f) corresponds to [/-statistics of order n = 1). Our investigation parallels 
the classical and well-developed theory of [/-statistics of independent random variables, however we would 
like to point out that in our context additional interest in this type of functionals of the process X stems 
from the fact that they capture 'average dependencies' between particles of the system. This will be seen 
from the form of the limit, which turns out to be more complicated than in the i.i.d. case. 

The organization of the paper is as follows. After introducing the basic notation and preliminary facts 
in Section [2] we describe the main results of the paper in Section [3j Next (Section U]) we restate the results 
in the special case of n = 1 (as proven in [1]) to serve as a starting point for the general case. Finally, in 
Section [5] we provide complete proofs for arbitrary n. We conclude with some remarks concerning the so 
called non-degenerate case (Section [6]). 

2 Preliminaries 
2.1 Notation 

For a branching system {Xj}^>q, we denote by \Xt\ the number of particles at time t, and by Xt{i) - the 
position of the i-th (in a certain ordering) particle at time t. We sometimes use E^. or P^. to denote the 
fact that we calculate the expectation for the system starting from a particle located at x. We use also E 
and P when this location is not relevant. 

By — )-'^ we denote the convergence in law. We use <, ~ to denote the situation when an equality or 
inequality holds with a constant c > 0, which is irrelevant for calculations. E.g. f{x) g[x) means that 
there exists a constant c > such that /(x) = cg{x). 



\Xt\ 





(2) 



il,i2v*n = l 

ife^ij if k^j 
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By X o y = XiUi we denote the standard scalar product of x,y £ W^, by || • || the corresponding 

Euchdean norm. By (8)" we denote the ?i-fold tensor product. 
We use also (/,//) := J^a f{x)fj,{dx). 

In our paper we will use the Feynman diagrams. A diagram 7 labeled by {1, 2, . . .} is a graph consisting 
of a set of edges not having common endpoints, and unpaired vertices A^. We will use r{'y) to denote 
the rank of the diagram i.e. the number of edges. For properties and more information we refer to \n\ 
Definition 1.35]. 

In the paper we will use the space 

p = -piW^) := |/ : M*^ h-^ M : / is continuous and 3„ such that ^ as ||3;|| +oo| , (3) 

We endow this space with the following norm 

\v = \\f\\v(m ■= \nix)f{x)\, (4) 



where n{x) := exp ■ We will also use 

to denote the space of continuous compactly supported functions. 

Given a function / G 7^(M'^) we will implicitly understand its derivatives (e.g. ^) in the space of 
tempered distributions (see e.g. |21[ p. 173]). 

2.2 Basic facts on the Galton- Watson process 

The number of particles {|-'^t|}t>o is the celebrated Galton- Watson process. We present basic properties 
of this process used in the paper. The main reference in this section is j^. In our case the expected total 
number of particles grows exponentially at the rate 

Xp := {2p - 1)A. (5) 

The process becomes extinct with probability (see [2l Theorem 1.5.1]) 

1-p 

Pe = • 

P 

We will denote the extinction and non-extinction events by Ext and Exf^ respectively. The process 
Vt := e~'^*'*|Xj| is a positive martingale. Therefore it converges (see also [21 Theorem 1.6.1]) 

Vt — ^ Voo, CL.s. as t — )• +00. (6) 
We have the following simple fact (we refer to [T] for the proof). 
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Fact 2.1. We have {Voo = 0} = Ext and conditioned on non- extinction V^o has the exponential distribution 
with parameter . We have E(VJx)) = 1 and Var(yoo) = 2p-i ■ IEe~^'^p*|Xf |^ is uniformly bounded, i.e. 
there exists C > such that for any t > we have E,e~'^^p^\Xt\^ < C. Moreover, all moments are finite, 
i.e. for any n G N and t > we have ElXjl" < +oo. 

We will denote the variable Voq conditioned on non-extinction by W. 
2.3 Basic facts on the Orstein-Uhlenbeck process 

We recall that the Ornstein-Uhlenbeck process is a time homogenous Markov process with the infinitesimal 
operator 

L := ^cr^A - i^ix o V. (7) 

The corresponding semigroup will be denoted by T. The density of the invariant measure of the Ornstein- 
Uhlenbeck process is given by 

V>ix) ■■= (-r^) exp(--^llxf) . (8) 



2.4 Basic facts concerning [/-statistics 

We will now briefly recall basic notation and facts concerning [/-statistics. A [/-statistic of degree n based 

on an Af- valued sample Xi, . . . , and a function / : X"^ — >• M, is a random variable of the form 

N 

ik^ij if k^j 

The function / is usually referred to as the kernel of the [/-statistic. Without loss of generality it can 
be assumed that / is symmetric i.e. invariant under permutation of its arguments. We refer the reader to 
[181 El foi' more information on [/-statistics of sequences of independent random variables. 

In our case we will consider [/-statistics based on the sequence of positions of particles from the 
branching system as defined by ([2]). We will be interested in weak convergence of properly normalised 
[/-statistics when i — t- oo. Similarly as in the classical theory, the asymptotic behaviour of [/-statistics 
depends heavily on the so called order of degeneracy of the kernel /, which we will briefly recall in Section 

A function / is called completely degenerate or canonical (with respect to some measure of reference 
9?, which in our case will be the stationary measure of the Ornstein-Uhlenbeck process) if 

/(Xi, . . .,Xn)^{dXk) = 0, 

for all xi, . . . , Xfc+i, . . . ,Xn ^ X. The complete degeneracy may be considered a centredness condi- 
tion, in the classical theory of [/-statistics canonical kernels are counterparts of centred random variables 
from the theory of sums of independent random variables. Their importance stems from the fact that each 
[/-statistic can be decomposed into a sum of canonical [/-statistics of different degrees, a fact known as 
the Hoeffding decomposition (see Section I5.1.2P . 



X 
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3 Main results 



This section is devoted to the presentation of our results. The proofs are deferred to Section [5j 

We start with the fohowing law of large numbers (throughout the article when dealing with [/-statistics 
of order n we wiU identify R^><_._><M^ with M"'^). 

n 

Theorem 3.1. Let {Xt}^^Q be the OU branching system starting from x ^ W^. Let us assume that 
f : M'"* I— 7- M is a bounded continuous function. Then, on the set of non- extinction Ext^ there is the 
convergence 

lim |Xir"C/f(/) = (/,(^®"> a.s. (9) 

Moreover, when f G V(W^'^), then the above convergence holds in probability. 

Having formulated the law of large numbers let us now pass to the corresponding CLTs. As already 
mentioned in the introduction, their form depends on the relation between Ap and fi, more specifically we 
distinguish three cases Ap < 2/x, Ap = 2/i and Ap > 2/i. We refer the reader to [T| for a detailed discussion 
of this phenomenon as well as its heuristic explanation and interpretation. Here we only stress that the 
situation for Ap > 2/i differs substantially from the remaining two cases, as we obtain convergence in 
probability and the limit is not Gaussian even for n = 1 (intuitively, this is caused by large branching 
intensity which lets local correlations between particles prevail over the ergodic properties of the Orstein- 
Uhlenbeck process). 

3.1 Slow branching case: Ap < 2/i 

Let Z he a Gaussian stochastic measure on M.'^'^^ with intensity /ii(dtdx) := ((5o(dt) + 2Xpe^p^dt) ip{x)dx 
defined according to \n\ Definition 7.17]. We denote the stochastic integral with respect to Z by / and 
the corresponding multiple stochastic integral by |17| Section 7.2]. We assume that Z is defined on 
the probability space (1^,7-", P). For / G 'P(]R"'^) we define (we recall that Tis the semigroup of the 
Orstein-Uhlenbeck process) 

H{f){si,Xi, S2,X2, ...,Sn, Xn) := i'^i=iTsi) f{xi,X2, Xn), Si G M+, Xi G R'^ . (10) 

It will be useful to treat this function as a function of n variables of type Zi := {si,Xi) G x R'^. For a 
Feynman diagram 7 labeled by {1, 2, . . . , n} 

^(/,7) :=^|A^| n J l^2idzj^k)\ H{f){ui,U2,...,Un)] , //2(dtdx) := 2ApeV(^(2;)dtdx, (11) 



Jj,k)eE. 



■7 



where Ui = Zj^k if (ji k) G and i = j or i = k and Uj = Zj if i G A^. Less formally, for each pair (j, k) 
we integrate over diagonal of coordinates j and k with respect to fi2- The function obtained in this way 
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is integrated using the multiple stochastic integral I\a.^\- We define 

Li(/):= J;L(/,7), (12) 

7 

where the sum spans over all Feynman diagrams labeled by {1, 2, . . . , n}. 

Fact 3.2. For any canonical f G V{W''^) we have ELi(/)^ < +oo. Moreover Li is a continuous function 
Li : Can i— t- L2(f^,-F, P), where Can = {/ G V{W^'^) : f is a canonical kernel^ and Can is endowed with 
the norm \\ ■ . 

We are now ready to formulate our main result for processes with the small branching rate. 

Theorem 3.3. Let {Xt}^^Q be the OU branching system starting from x S W^. Let us assume that 
f S P(M"'^) is a canonical kernel and Xp < 2fi. Then conditionally on the set of non- extinction Exf^ there 
is the convergence 

(^e-^^*l^d,^^^=p^,^) -^'{W,G„L,{f)), (13) 
where where Gi ~ A/'(0, l/(2p — 1)) and W, Gi, Li{f) are independent random variables. 

3.2 Critical branching case: Xp = 2/i 

Consider the space C := L2(M^,$(dx)), where 

\d(p{x) 



dxi 



$(x) := 

1=1 

and a centred Gaussian process {Gf)f^c defined on some probability space (0, J^, P), with the covariance 
structure given by 

c„v,G,.G,).^E(/.S){/.S)- 

^1=1 ' ' 

We will identify the process with a map 1: C ^ L2(0,J-', P), such that /(/) = Gj. One can easily 
check that / is a bounded linear operator. 

To formulate the central limit theorem in this case we will need the following 

Fact 3.4. For every n > 1 there exists a unique bounded linear operator L2 : Can 1— )• ^2(^2, P), 
where Can = | / G P(M"'^) : f is a canonical kernel^ endowed with the norm \\ ■ \\^, such that for every 
/i, ...,/„: ^ M satisfying {fi, if) = and fi G V{W^), i = 1, . . . , n, 

^2(/i®---® /„)=/(/!)••• /(/„). (15) 

The above lemma will be proved together with the following theorem, which describes the asymptotic 
behaviour of ?7-statistics in the critical case. 
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Theorem 3.5. Let {Xj}j>q be the OU branching system starting from x G . Let us assume that 
f S V{W^'^) is a canonical kernel and Xp = 2fj,. Then conditionally on the set of non- extinction Exf^ there 
is the convergence 

where Gi ~ A/'(0, l/{2p — 1)) and W, G, L2{f) are independent random variables. 
3.3 Fast branching case: Ap > 2/i 

In order to describe the limit we introduce an M'^-valued process 

Ht := e(-^''+^)* Yl ^*(^)' * ^ 0- (16) 

i=l 



pd 



The following two facts have been proved in [T]. 

Fact 3.6. H is a martingale with respect to the filtration of the OU-branching system starting from x G 
Moreover for Xp > 2^ we have sup^Ei^^ < +oo, therefore there exists := limf^^ca Ht (a.s. limit) and 
Hoa G When the OU branching system starts from then martingales Vt and Ht are orthogonal. 

It is worthwhile to note that the distribution of Hqo depends on the starting conditions. 

Fact 3.7. Let {X(}^>q and {Xt}t>Q be two OU branching processes, the first one starting from and the 
second one from x. Let us denote the limit of corresponding martingales by Hoo, H^o respectively. Then 

where Voo is the limit given by ([6]) for the system X. 

Hao is M'^- valued, we denote its coordinates by W^. Let / S V{MP''^). We define 

L-^(f)-= V / — ...J^in 

where we adopted convention that xj^i is the l-th. coordinate of the j-th variable. By L^[f) we will denote 
^sif) conditioned on Exp^. 

Theorem 3.8. Let {Xt}^^^ be the OU branching system starting from x G M.'^. Let us assume that 
f S 7^(R"''^) is a canonical kernel and Xp > 2fj,. Then conditionally on the set of non- extinction Exf^ there 
is the convergence 

|^e-V|Xi|,^iL_0!:°£,e-n(A,-M)*[/n(/)j _,d (^,Gi,L3(/)), (17) 

where Gi ~ A^(0, l/(2p — 1)) and {W,L3{f)),Gi are independent. Moreover 

(e-^-*|Xt|,e-"(^^-'^)*f/,"(/)) -^'^ (v^,Ls{f)^ , in probability. 
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3.4 Remarks on the CLT for [/-statistics of i.i.d. random variables 

For comparison purposes we will now briefly recall known results on the central limit theorem for U- 
statistics of independent random variables, [/-statistics were introduced in the 1940's in the context 
of unbiased estimation by Halmos |14| and Hoeffding who obtained the central limit theorem for non- 
degenerate (degenerate of order 0) kernels |16) . The full description of the central limit theorem was 
obtained in \10\ (see also the article [T3] where the CLT is proven for a related class of F-statistics). 
Similarly as in our case, the asymptotic behaviour of [/-statistics based on a function /: Af" — )■ M and an 
i.i.d. A'-valued sequence Xi^X2, ... is governed by the order of degeneracy of the function / (see Section 
I5.1.2P with respect to the law of Xi (call it P). The case of general / can be reduced to the canonical one, 
for which one has the weak convergence 

TV 

Ar-n/2 ^ /(X,,,...,X,^)^ J„(/), 

jl,j2v«n = l 

jfe^ij if kytj 

where J„ is the n-fold stochastic integral with respect to the so-called isonormal process on Af, i.e. the 
stochastic Gaussian measure with intensity P. 

Let us note that in the i.i.d. case the limiting distribution is simpler than for [/-statistics of the 
OU branching processes. For small branching rate however, the behaviour of [/-statistics in our case 
still resembles the classical one as it is a sum of multiple stochastic integrals of different orders. In 
the remaining two cases the behaviour differs substantially. This can be seen as a result of the lack of 
independence. Although asymptotically the particles' positions become less and less dependent, in short 
time scale offspring of the same particle stay close one to another. 

Let us finally mention some results for [/-statistics in dependent situations, which have been obtained 
in the last years. In [6] the authors analysed the behaviour of [/-statistics of stationary absolutely regular 
sequences and obtained the CLT in the non-degenerate case (with Gaussian limit). In [5] the authors 
considered a and ip mixing sequences and obtained a general CLT for canonical kernels. Interesting 
results for long-range dependent sequences have been also obtained in |8]. A more recent interesting work 
is |19) . where the authors consider [/-statistics of interacting particle systems. 

4 The case of n = 1 

In the special case of n = 1 the results presented in the previous section were proven in [1] . Although this 
case obviously follows immediately form the results for general n it is actually a starting point in the proof 
of the general result (similarly as in the case of [/-statistics of i.i.d. random variables). Therefore, for 
the reader's convenience, we will now restate this case in a simpler language of [1], not involving multiple 
stochastic integrals. 

We will start with the law of large numbers 
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Theorem 4.1. Let {Xt}^^Q be the OU branching system starting from x G M.'^. Let us assume that 
f G V{R'^). Then 

lim e""^*"* {Xt, f) = (/, (f) Voo in probability, 

t— S- + CXD 

or equivalently on the set of non- extinction, Ext^, we have 



lim \Xt\ ^ {Xt, f) = {f, If) in probability. (18) 

t— >-+oo 



Moreover, if f is bounded then the almost sure convergence holds. 

4.1 Small branching rate: Xp < 2/i 

We denote f{x) := f{x) — {f,<f) and 



a 



f 



^, p) +2Xp [ °° (v, (e^^^'^^Tsf)) ds. (19) 







Let us also recall ([6]) and that W is Voo conditioned on Ext^. In this case, the behaviour of X is given by 
the following 

Theorem 4.2. Let {Xt}^^^ be the OU branching system starting from x G M'^. Let us assume that 
Xp < 2fi and f G V{W^). Then < +oo and conditionally on the set of non- extinction Exf^ there is the 
convergence 

where G\ ~ A/'(0, l/(2p — 1)), ~ A/'(0, CTj) and W, Gi, G2 are independent random variables. 

4.2 Critical branching rate: Xp = 2/i 

We denote 



1=1 



Note that the same symbol has already been used to denote the asymptotic variance in the small 
branching case. However, since these cases will always be treated separately, this should not lead to 
ambiguity. 

Theorem 4.3. Let {Xj}j>q be the OU branching system starting from x G M.'^. Let us assume that 
Xp = 2fi and f G V{W^). Then cJj < +00 and conditionally on the set of non- extinction Exf^ there is the 
convergence 

where Gi ~ AA(0, l/(2p — 1)), G2 ~ AA(0, 0"^) and W, Gi, G2 are independent random variables. 
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4.3 Fast branching rate: \p > 2/i 

In the following theorem we use the notation introduced in Section 13.31 



Theorem 4.4. Let {Xt}^^^ be the OU branching system starting from x G W^. Let us assume that Xp > 2/i 
and f G 'P(R'^). Then conditionally on the set of non- extinction Exf^ there is the convergence 

1^^ l''^'' ^ ' exp((A-/i)t) WG,(grad/,^)oJ), (21) 

where G ~ AA(0, l/(2p — 1)), (VF, J), G are independent and J is Hoo conditioned on Ext^. Moreover 

e-V|X,|, iXtJ)- \Xt\{f ,ip) \ ^ ^^^^ ^g^^^^^ ^ ^^^^ probability. 
exp((A-/x)t) y 

5 Proofs 

We will now pass to the proofs of the results announced in Section [3l Their general structure is similar as 
in the case of [/-statistics of independent random variables, i.e. all the theorems will be proved first for 
linear combinations of tensor products and then via suitable approximations extended to proper function 
spaces. 

In the next section we will recall some additional (standard) facts concerning the Orstein-Uhlenbeck 
process and [/-statitistics. Next, in Section [^T^ we will develop general tools needed for the approximation, 
which will be the most technical part of the proof. Section 15.31 will be devoted to rather short proofs of 
the main results. 

From now on we will often work conditionally on the set of non-extinction Ext^, which will not be 
explicitly mentioned in the proofs (however should be clear from the context). 

5.1 Auxihary facts and notation 
5.1.1 The Orstein-Uhlenbeck process 

The semigroup of the Ornstein-Uhlenbeck process can be represented by 

Ttf{x) = {gt*f){xt), xt:=e-^*x, (22) 

where 

9ti^) = {-^lY' exp I , := a\l - e'^A'*). 

Let us recall ([8]). We denote ou{t) := i/l — e"^^* and let G ip. The semigroup Thas the following useful 
representations 

Ttf{x)= [ f{xt-y)gt{y)dy= [ / (xe"^* + on(t)y) (^(y)d2/ = E/(xe-^* + on(t)G). (23) 



We also denote 



r/ := e^'Ts. 
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5.1.2 U- and ^/-statistics 



We will now briefly recall one of the standard tools of the theory of fZ-statistics, which we will use in the 
sequel, namely the Hoeffding decomposition. 

Let us introduce for / C {1, . . . ,n} the Hoeffding projection of / corresponding to / as the function 
11//: RI^I'^ — )• R, given by the formula 

Uifiixi),^j) = [ (ll^idy,)ll{5,^-ip)idyi))fiyi,...,yn). (24) 

Once can easily see that for |/| > 1, 11// is a canonical kernel. Moreover 110/ = J^^a /{xi, ■ ■ ■ , Xn) ITiLi f{dxi). 

Note that if / is symmetric (i.e. invariant with respect to permutations of arguments), 11// depends 
only on the cardinality of /. In this case we speak about the k-th Hoeffding projection {k = 0, . . . ,n), 
given by 

nfc/(xi, ...,Xk) = {{6,, - if) (S> ■ ■ ■ 'S) (<5,, -ip)<S> /) . 

A symmetric kernel in n variables is called degenerate of order k — 1 {1 < k < n) iS k = min{i > 
0: Hif ^ 0}. The order of degeneracy is responsible for the normalisation and the form of the limit in 
the central limit theorem for [/-statistics, e.g. if the kernel is non-degenerate, i.e. Hi/ ^ 0, then the 
corresponding [/-statistic of an i.i.d. sequence behaves like a sum of independent random variables and 
converges to a Gaussian limit. The same phenomenon will be present also in our situation (see Section [6]). 

In the particular case k = n the definition of the Hoeffding projection reads as 

n„(/) := m=i{s.. - ^)) J) ■ 

One easily checks that 

f{xi,...,Xn) = ^ n//((xi)i6/), 
7C{l,...,n} 

which gives us the aforementioned Hoeffding decomposition of [/-statistics 

/C{l,...,n} ^' ' 

which in the case of symmetric kernels simplifies to 

-r(/)^|:(;)M.?(n./), 

where we use the convention U^{a) = a for any constant a. 

For technical reasons we will also consider the notion of a ^/-statistic which is closely related to U- 
statistics, and is defined as 

yru)-= E f{Xt{ii),xt{i2)....,xt{in)). (25) 

11 



The corresponding Hoeffding decomposition is 

vr{f)= i^ir-'"v;'"(n//), (26) 

/C{l,...,n} 

where again we set Vf'{a) = a for any constant a. 

In the proof of our results we will use a standard observation that a ?7-statistic can be written as 
a sum of y-statistics. More precisely, let J' be the collection of partitions of {1, . . . ,n} i.e. of all sets 
J = {Ji,...,Jfc}, where Jj's are nonempty, pairwise disjoint and {Jj^Ji = {l,...,n}. For J as above 
let fj be a function of |J| variables xi, . . . , x^j^, obtained by substituting Xi for all the arguments of / 
corresponding to the set Jj, e.g. for n = 3 and J = {{1,2}, {3}}, fj{xi,X2) = f{xi,xi,X2)- An easy 
application of the inclusion-exlusion formula yields that 

where aj are some integers depending only on the partition J. Moreover one can easily check that if 
J = {{1} , . . . , {?^}}, then aj = 1, whereas if J consists of sets with at most two elements then aj = (—1)'^ 
where k is the number of two-element sets in J. Let us also note that partitions consisting only of one- 
and two-element sets can be in a natural way identified with the Feynman diagrams (defined in Section 

EH). 

5.2 Approximation 

5.2.1 Approximation of functions 

First we will show that any function in V{W^'^) can be approximated by tensor functions. For a subset A 
of a linear space by span{A) we denote the set of finite linear combinations of elements of A. 

Lemma 5.1. Let A := {&l^ifi '■ fi bounded continuous} and f € V{W^'^) be a canonical kernel then there 
exists a sequence {fk} C span{A) such that each fk is canonical and 

ll/fc - /Hp 0, ask^+oo. 

Proof. First we prove that span{A) is dense in T'(M"'^). Let us notice that given a function / G 'P(M"''^) it 
suffices to approximate it on some box [— M, M]'^, M > 0. The box is a compact set and an approximation 
exists due to the Stone- Weierstrass theorem. 

Now, let / G 'P(M"'^). We may find a sequence {hk} C span{A) such that /i^ — )• / in "P. Let us recall 
the Hoeffding projection (j24p and denote / = {l,2,...,n}. Now direct calculation (using the exponential 
integrability of Gaussian variables) reveals that the sequence fk ■= H/Zifc fulfils the conditions of the 
lemma. □ 
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Let / G •P(M"'^) and / C {1, 2, . . . , n} with / = k. We define 

fl{xi,X2,. . . ,Xn) ■■= / f{zi,Z2,...,Zn)Y\9liXie~ 



yi)dyi, 



(28) 



where Zj = yj if i G /, = Xj otherwise and gi is given by (p2|) . We have 

Lemma 5.2. Lei / € 'P(]R"''^) and / € N. Then for any I C {1, 2, . . . , n} i/ie function fj is smooth with 
respect to coordinates in I. For any multi-index A = C {1, . . . ,ndy such that {[ij/d]: j = 

1, ... ,1} C I we have 

11-5^11. <c||/||„ 

where C > depends only on a, fi,d,l,n. Moreover, when f is canonical, so is //. 



/ f{zi,Z2, Zn)- — TTsila^je"'" - yi)dyi. 

jRkd dXA ^J- 



Proof. Let us fix some / and A. Using (j28|) we get 



OXA 



Therefore by the properties of Gaussian density gi and easy calculations we arrive at 

aiAi 



\K\ < 



■p / n(zi,Z2, ■■■,z„ 



dxA J 



dyi < CAn{xi,X2, ■■■,Xn) \ 



for some constant Ca- 

To conclude it is enough to take the maximum over all admissible pairs /, A. 

Let us now assume that / is canonical. We would like to check that for any j G {1, 2, . . . , n} we have 

/ f{xi,X2, . . . ,Xn)(p{Xj)dXj = 0. 
There are two cases, the first when j ^ /. Then we have 

/ f{xi,X2,...,Xn)f{xj)dxj= / f{zi,Z2,. ■ ■ ,Zn)\\gi{xie~^ - yi)dyiLp{zj)dzj 

= (/ fizi,Z2,...,Zn)(pizj)dzj]T]gi{xie~'' -yi)dyi = 0. 

The second case is when j £ I. Then 

f{xi,X2, . . . ,Xn)^{Xj)dXj 

= f{zi,Z2,...,Zn) TT giixiC"'^ - yi)dyi I giixje^f" - yi)ip{xj)dxj] dyj 

Jr"" .g^^^.^ KJr" j 

= f{zi,Z2,...,Zn) TT giixie''' - yi)dyiip{yj)dyj, 

where the second equality holds by the fact that (p is the invariant measure of the Ornstein-Uhlenbeck 
process. Now the proof reduces to the first case. □ 
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We will also need the following simple identity. We consider {xj}-^-^ ^ n > {^i}i=i 2 n- induction 



one easily checks (we slightly abuse the notation here, e.g. ^ denotes the derivative in direction 
and the integral over the segment [a,b] C M*^) 

Lemma 5.3. Let f be a smooth function, then 

+ ei(xi - Xi),X2 + e2{x2 - ^2), . . . + en{Xn - Xn)) = 

(ei,e2,...,e„)e{o,ir 



Xi Xn 



gn 



dyidy2 ...By, 



-fiyi,y2, ■ ■ . , yn)dy„dy„_i . . .dyi. 



5.2.2 Approximation of J7- statistics 

Bookkeeping of trees So far we have shown that one can approximate functions in V by linear com- 
binations of tensor products. Our next goal is to show that two functions which are close in V generate 
?7-statistics which (after proper normalization, specific for each regime) are close in distribution. To this 
end we will use the "bookkeeping of trees" technique (see e.g. [3] or [101 Section 2]), which via some 
combinatorics and introduction of auxiliary branching processes will allow us to pass from equations on 
the Laplace transform in the case of n = 1 to estimates of moments of ^/-statistics and consequently 
[/-statistics. 

We recall 1^. Let /i, /2, ...,/„ £ C^iR"^) and fi > 0. We would like to calculate 

/ n \ 

(29) 



w(®r=i/.) = EJn^^*'/*: 



\i=l 



Let A C {1, 2, . . . , n}, slightly abusing notation we denote 
WA{x,t,a) := Ex ^exp j-^Oj , 



WA{x,t) := - — w{x,t,a) 



a=0 



Note that this differentiation is valid by Fact 12. II and properties of the Laplace transform (e.g. |12^ Chapter 
XIIL2]). By the calculations from Section 4.2. in [1] we know that 

t(;A(x,t,a) = 7Ie"^'eA"«/.(x) + A / Tts [pwl{-, s,a) - wa{-, s,a) + {I - p)] ds. 

Jo 

It is easy to check that 

(-l)l^lu;A(x,t). (30) 



VA 



{x,t) :=E,. (n(^t./i> J 

ViGA / 



Assume that |A| > 0. We denote by -Pi (A) aU pairs (Ai, A2) such that Ai U A2 = A and Ai n A2 = 0, and 
by P2{^) C -Pl(A) pairs with an additional restriction that Ai 7^ and A2 7^ 0. We easily check that 



daA 



WA 



\ieA / 



+ A / Tt 



p 



E 



(Ai,A2)ePi(A) 



glAil c}|A2| ^|A| 

wa{-, s, a)- wa{-, s, a) - - — wa{-, s, a) 

OaAi aOAa OUA 



{x)ds. 
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We evaluate it at a = 0, (let us notice that g^^' wk{x, s, 0) 
sides by (—1)1^1 and use the definition of v\{x,t) to get 



§^'WAi{x,s,0) =WAi{x,s)), multiply both 



{^,t) = rtlllh] (x) + A [Wt^ 
= rt(llfi] (^) + ^ f'T't- 

ViGA / 



P X] VA-,{-,s)vA2{-,s) - va{-,s) 

(Ai,A2)GPi(A) 



(x)ds 



P X] VA^{-,s)vA2{-,s) + {2p - 1)va{-,s) 
(Ai,A2)GP2(A) 



(x)ds. 



This can be easily transformed to (recall that Ts'^ f{x) = e'^p^Tsf{x)) 



VA{x,t) = Tt 



VA,{-,S)VA2{-,S) 

{Ai,A2)GP2(A) 



(x)ds. 



(31) 



The last formula is much easier to handle if written in terms of auxiliary branching processes. Firstly, we 
introduce the following notation. For n G N \ {0} we denote by T„ the set of rooted trees described below. 
The root has a single offspring. All inner vertices (we exclude the root and the leaves) have exactly two 
offspring. For r G T.„, by 1{t) we denote the set of its leaves. Each leaf I £ 1{t) is assigned a label, denoted 
by lab(/), which is a non-empty subset of {1, 2, . . . , n}, the labels fulfil two conditions: 



IJ lab(0 = {l,2,...,n}, Vz,,i,6;(,) (/i / /2 
m{T) 



lab(/i)nlab(/2) 



In other words, the labels form a partition of {1,2, . . . ,n}. For a given r G T„ let f(r) denote the set of 
inner vertices (we exclude the root and the leaves), clearly |«(t)| = \1{t)\ — 1 (as usual | • | denotes the 
cardinality). Let us identify the vertices of r with {0, 1, 2, . . . , |r| — 1} in such way, that for any vertex i 
its parent, denoted by p(i), is smaller. Obviously, this implies that is the root and that the inner vertices 
have numbers in set {1, 2, . . . , }. We denote also 



s(r) := {I E 1{t) : | lab(/)| = 1} , m(T) := {I £ 1{t) : \ lab(/)| > 1} , 



(32) 



leaves with single and multiple labels respectively. 

Let r € T„, we consider an Ornstein-Uhlenbeck branching walk on r as follows (in this part we ignore 
the labels). Let us fix t G and {ti}^^^(^^y The initial particle is placed at time at location x, it evolves 
up to the time t — ti and splits into two offspring, the first one is associated with the left branch of vertex 
1 in tree r and the second one with the right branch. Further each of them evolves until time t — ti, where 
i is the first vertex in the corresponding subtree, when it splits and so on. At time t the particles are 
stopped and their positions are denoted by (the number of particles at the end is equal to the 

number of leaves). The construction makes sense provided that ti < t and ti < tp(^i-^ for all i G ^(t). We 
define 



OC/(»r=i/^,T,t,{tO,g,(,),x) :=E m/,(y,(,)) , 



(33) 



\a=l 
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where j{a) = I £ 1{t) is the unique leaf such that a £ lab(/). We also define 



SiT,t,x) := (pA)l'MleV 



n 

I iGfi(r) 



(34) 



Fact 5.4. Let A„ = {1, 2, . . . , n}. T/ie following identity holds 



VAn{x,t) = ^ S'(T,t,x). 



(35) 



Proof. The claim is a consequence of the identity v\{x,t) = ^^^-^^(A) S{T,t,x), where A C {1,2, . . . , n} 
and T(A) is the set of trees, as Tr|^|, with the exception that the labels are in the set A. This identity 
in turn will follow by induction with respect to the cardinality of A. For A = {i} equation ()3ip reads as 
VA{x,t) = e'^p^Ttfi{x). The space T(A) contains only one tree, denoted by r^, consisting of the root and a 
single leaf labelled by i. We have i(r) = and obviously S{Ts,t,x) = e^p^Ttfi{x). Let now |A| = /c > 1. 
Similarly as before the first term of (j3ip corresponds to r^. By induction the second term can be written 
as 



Jo 



E E E S{n,t^,■)S{T2,t^,■) 

(Ai,A2)GP2(A) riGT(Ai) T2GT{A2) 



{x)dti. 



We have 



(Ai,A2)GP2(A) riGT(Ai) r2eT(A2) 



e-^^'^p{t-ti,x,y) 



dti 



dti. 



I iei(ri) 



y-iei(T2) 



dydi 



1, 



where p is the transition density of the Ornstein-Uhlenbeck process. Now we create a new tree r by setting 
Ti and T2 to be descendants of the vertex born at time t — ti. We keep labels and split times unchanged 
and assume that the first particles of OU processes on ri and T2 are put at x. Thus by the Markov 
property of the Ornstein-Uhlenbeck process we can identify the branching random walk on ti and T2 with 
the branching random walk on r. It is also easy to check that the described correspondence is a bijection 
from set of pairs (ti, r2) (as in the sum above) to T„ \ {r^} and therefore the expression above is equal to 
the sum ET„\{r4 S{T,t,x). □ 

The calculations will be more tractable when we derive an explicit formula for {^j}jg/(T-)- Let us recall 
the notation introduced in (|23|) and consider a family of independent random variables {Gi}-^^, such that 
Gi ^ ip for i 7^ and Gq ~ Recall also that ou{t) = \/l — e"^'^*. The following fact follows easily from 
the construction of the branching walk on r and ()23p . 
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Fact 5.5. Let be positions of particles at time t of the Ornstein-Uhlenheck process on tree r with 

labels {ti}i(zi(j) ■ We have 

where for any i G Z(r) we put 



Zi:= ^ ou{tp(i) - ti)Gie + ou{tp(^i))Gi, 

where P{i) := {predecessors of i}, by convention we set to = t and ou{tp(^Q) — to) = 1- 
We are now ready to prove an extended version of Fact [57 



Fact 5.6. Let {X(}^>q be the OU branching system starting from x G M'^ and f G V{W^''^) then 

E J2 fiXtiii),Xtii2),...,Xt{in))= S{T,t,x), (36) 
where in ([M]) we extend the definition of OU in (j33p by putting 

ou{f, T, t, {t , x) := E/(y,-(i), y,(2), . . . , y,(„)). 

Moreover all the quantities above are finite. 

Proof. (Sketch) Using Fact I5.4| Fact 12.11 (|29p and the Lebesgue monotone convergence theorem together 
with hnearity one may prove that (f36|) is vahd for / = C, C > 0. Using standard tricks we may drop the 
positivity assumption in (p9|) and ()35p . Therefore, by the Stone- Weierstrass theorem and the dominated 
Lebesgue theorem, dMI) is vahd for any / £ C^{W"^). Let now / E V{W"^), / > 0. We notice that for 
any r € T„ the expression OU{f,T,t,{ti}-^-^^^ ,x) is finite, which fohows easily from Fact 15.51 Further, 
one can find a sequence {fk} such that fk € C^{M.'^), //c > and fk f (pointwise). Appeahng to the 
monotone Lebesgue theorem yields that (j36p still holds (and is finite). To conclude, once more we remove 
the positivity condition. □ 

As a simple corollary we obtain 

Corollary 5.7. Let {Xt}^^^ be the OU branching system, then for any n > 1 there exists Cn such that 

Proof. We apply the above fact with / = 1. Using definition (I34p and the inequality \i{T)\ < n — 1 
for r € T„, it is easy to check that for any t G T„ we have S{T,t,x) < Ct-c"'^*'*, for a certain constant 
depending only on r. □ 

Let us recall notation of (1321). We have 
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Fact 5.8. For any n G N there exists C,c > 0, such that for any r G T„ any {ti}jgj(,-), split times as 
above, and any canonical f G V{W^'^) we have 



OC/(/,r,t,{ti}.g.(^),x) < ||/||pCexp{c||x||}exp(^-/i ^ tp(^i)j. (37) 

i£s{T) 

Proof. Let k < n. Without loss of generality we may assume that / := {1, 2, . . . ,k} are single numbers 
(i.e. j{i) G s{t)) and {k + 1, . . . ,n} are multiple ones. Let us also assume for a moment that for i G 
{1, 2, . . . ,k} we have tp[j[i)) > 1- Let Zj and Gi be as in Fact 15.51 We have E/(yj(i), 1^(2), • • • , ^(n.)) = 
E/(Zj(i), Zj(2), . . . , Zj(n))- For i < we define 

Zi := ou{tp^i) - ti)Gie-^^'^-^^ + ou{tp^i) - 1)G,. 
Moreover, let / := // be given by (j28p . By the semigroup property of Twe have 

^f{Zj(l), Zj{2)-, ■■■ , Zji^n)) = lE/(Zj(l), . . . , Zj(fc), . . . , Zj(ri)) =: A. 

By the fact that / is smooth with respect to coordinates from / and canonical and Lemma 15.31 we obtain 
^ = IE ^ (-l)^'=i^'/(Zj(i) + ei(Gj(i) - ^j(i)), • • • , -^j(fc) + efc(G'j(fc) - -^j(fc)), ^j(fc+i), • • • , ^i(n)) 

(ei,...,efc)e{0,l}'" 

= 1E/ ... / — — /(yi,...,yfc,Z^(fc+i),...,Zj(„))dyfc...dyi. 

JGj(i) ■'Gj(fc) (^yi • • • <^t/fc 

From now on we restrict to the case d = \. The proof for general d proceeds along the same lines 
but it is notationally more cumbersome. Using Lemma 15.21 and applying the Schwarz inequality multiple 
times we have 



1^1 <C^II/llpE 



rZ. 



n / exp{|yi|}dyi j JJ exp{|Zj(i)|} 

<ll/llpnil/ exp{|yi|}dyi||2. JJ || exp II2, 



1=1 i=k+l 



Note that by the definition of Zi we have Zi — Gi = i7je~^^*p(^)~^^ + [ou{tpi^i) — 1) — l) Gj, where Hi 
is independent of Gi and i^j ~ J\f{xi,af) with fjj < a / y/2fl and ||xj|| < ||x||. Thus Zi — Gj is a Gaussian 
variable with the mean bounded by C||xj||e~^*J'(') and the standard deviation of order e~^^p^^K In particular 
\\Zi - Gi\\k < Gfcexp(Gfc||x|| - /utp(i)). Since 



/ '"exp{|yi|}d2/, < (el^^wl +el^^wl)|Zi 



the proof can be concluded by yet another application of the Schwarz inequality and standard facts on 
exponential integrability of Gaussian variables. 
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Finally, if some i's do not fulfil tp(j) > 1 we repeat the above proof with s(r) replaced by the set s' of 
indices for which additionally > 1. In this way we obtain (j37p with X^igs' ^p(i)- our setting 

hence ([37]) still holds (with a worse constant C). □ 

Bounds for U- and ^-statistics Using Fact 15.81 we will now be able to estimate moments of V- and 
{/-statistics of the branching particle system. The inequalities will be stated separately for each regime 
discussed in Section [3j 

We will first develop L2 bounds for [/-statistics with deterministic normalization. Since in the slow 
and critical branching case the normalization in our theorems is random, related to the size of the process, 
we will later transform these bounds to Lo-bounds for [/-statistics with a random normalization. 

Fact 5.9 (Small branching rate). Let {Xt}j.^Q be the OU branching particle system with Xp < 2fi. There 
exist C,c > such that for any canonical kernel / G "P (M"'^) we have 

2 



Proof. We need to estimate E,.e-"^^* ElS,...i2„=i fiMh), • • • , Xt{in))f{Xt{i2), • • • , Xt{i2n))- Obviously 
the function / (X" / is canonical. Moreover, it is easy to check, that ll/^/H-p < ||/||^- 

By Fact 15.61 it suffices to show that for each r G there exist C, c > such that for any i > we 
have e-"V5(r,t,x) < Cexp {c||2;||} 

Let us fix T G and denote by Pi{t) and P2iT) the sets of inner vertices of r with respectively one 
and two children in s{t). Set also Psir) := ^(t) \ (-Pi(t) U P2 (''")) • 

By the definition of S{t, t, x), Fact 15.81 and the assumption Xp < 2fi, we get 

e-^^pts(T,t,x) 

<Cie-("-i)V I JJ Tdt, I I n e'^*' J Of/ (/,r,t, X 
<C2||/||^exp(c||x||)e-("-i)^''* I [] f I U e^"'' ^U<t,,,, e-'^V.) 

|2 



<C2\\f\\i, exp{c\\x\\) exp ^{-{n - 1) + |Pi(t)|/2 + |P3(r)|)ApiJ 



X 

iGPilT)""" iGP2(T)"'" iGP3(r) 

<CM\v exp(c||x||)exp ((-(n - 1) + |Pi(t)|/2 + \Pz{r)\)Xpt 
To end the proof it is thus sufficient to show that |Pi(r)| -|- 2|P3(r)| < 2n — 2. 
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Note that |Pi(r)| + 2|P2(r)| = |s(r)| and 1^(^)1 = NWI = \Kr)\ - 1- Thus |Pi(t)| + 2|P3(t)| = 

2|Z(r)| — 2 — \s{t)\ = |/(r)| + |?TT.(r)| — 2 < 2n — 2 (recah that m(r) denotes the set of leaves with muhiple 
labels). This ends the proof. □ 

We will now pass to an analogous estimate in the critical case. 

Fact 5.10 (Critical branching rate). Let {Xt}^^Q be the OU branching particle system with \p = 2/i. 
There exist C, c > such that for any canonical kernel f G V{W''^) we have 

2 



Ex (^t-«/2g-(n/2)Aptyn^_^)^ < Cexp{c||x| 



Proof. We will use similar ideas as in the proof of Fact 15.91 as well as the notation introduced therein. 
Consider any r G T2n- By the definition of S{T,t,x), Fact 15.81 and the assumption Xp = 2//, we obtain 



t-"e-"^p*5(r,t,x) 
<C7i 11/11^ exp{c||x||}t-"e-('"-^)^''* 



<C2\\ffv exp(c||x||)t-" exp ((-(n - 1) + |Pi(t)|/2 + \P^{T)\)\pt ^ 
<C2||/||^exp(c||x||), 

where we used the fact that |P2(''")| < n and the estimate |Pi(t)| +2|P3(t)| < 2n — 2 obtained in the proof 
of Fact \5M 

□ 

Fact 5.11 (Fast branching rate). Let {Xt}^^^ be the OU branching particle system with Xp > 2/i. There 
exist C, c > such that for any canonical kernel f G 'P(M"'^) we have 



fe-"(^--'^)V,"(/)) < C7exp{c||x||} 
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Proof. As in the previous cases, consider any r G T2n- We have 
<C2||/||^exp(c||x||)e-2"(^^-'^)*+^^* 

<C^3||/||p exp(c||x||) exp ( - 2n(Ap - + Apt + |Pi(r)|(Ap - + |P2(r)|(Ap - 2^)t + jPsWIApt) . 

Thus it is enough to prove that 

Ap + |Pi|(Ap + \P2\{\p - 2fi) + iPslAp < 2n(Ap - (38) 

where for simphcity we write Pi instead of Pi{T) (in the rest of the proof we wiU use the same convention 
with other characteristics of r). Using the equahty |s| = |Pi| + 2|P2|) we may rewrite p8|) as 

Ap + |s|(Ap - /i) - |-P2|Ap + l-PslAp < 2n(Ap - fi), 

so by the inequahties 2n > |s| and Ap — /i > Ap/2 it is enough to prove that 

2 + |s| - 2IP2I + 2IP3I < 2?i. 

But IP3I = \i\ — IP2I — l-Pil and so 

2 + |s| - 2IP2I + 2IP3I = 2 + |s| + 2\i\ - 4IP2I - 2|Pi| = 2 - |s| + 2\i\ = -\s\ + 2\l\ = 2\m\ + |s| < 2n, 
which ends the proof. □ 

Random normalization Using the facts obtained above we wiU now prove estimates for [/-statistics 
normahzed by a proper power of \Xt\, which wih be relevant in the proofs of Theorems 13.31 and 13.51 Since 
asymptoticaUy \Xt\ exp(— Apt) behaves (conditionahy on Ext^) as an exponential random variable W and 
EVl^~^ does not exist, we will have to introduce a truncation, cutting out the set where \Xt\ is small. 

Corollary 5.12. Let {Xt}^^Q be the OU branching system with Ap < 2/i. There exist constants C,c > 
such that for any canonical f G P(]R"'^) and r G (0, 1) we have 

(w|Xi|-"/Vr(/)|l{|x.|>re^.'}) < Cexp{c||x||}r-/2||/(2n-)||p. 
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Proof. Let J' be the collection of partitions of {1, . . . ,n} i.e. of all sets J = { Ji, . . . , J^}, where Jj's are 
nonempty, pairwise disjoint and Jj = {1, . . . ,n}. Using (j27p and notation introduced there we have 

\x,\--/^ur{f) = ajvy^ifj), (39) 

Jej 

where aj are some integers depending only on the partition J. Since the cardinality of J depends only 
on n, it is enough to show that for each J £ J and some constants C, c > we have 

(|Xr"/'|^J"(/j)|l{|x.|>re^.'}) < C7exp{c||x||}r-"/2||/(2n.)||p. 

Let us thus consider J = { Ji, . . . , J^} and let us assume that among the sets Jj there are exactly I sets 
of cardinality 1, say Ji, . . . ,Ji. We would like to use Fact 15.91 To this end we have to express V^{fj) as a 
sum of y-statistics with canonical kernels. This can be easily done by means of Hoeffding's decomposition 
(j26|) . Since fj is already degenerate with respect to variables rri, . . . we get 

\Xt\-^'^v^\fj)= l^tl'"-"/'^J'''(n/c/j), (40) 

/C{Z+l,...,fc} 

where I*^ := {1, . . . , A:} \ /. Let us notice that n > 2A; — so n — A; > A; — / > |/|, which gives 

n-2|/| > A;- |/| = (41) 

Thus we have 

(|X,|l^l-"/2|yl^^l(n,c/,)|l^l^^l^^^.,.^) < (e-^^*("-2|/|)/2|^|/1(n,./,)|) 

for / 7^ {1, . . . , A;}, where in the third inequality we used Fact 15.91 One can check that for any n>2 there 
exists C > such that for any I,J C {1,2, ...,n} we have ||n//||-p < C||/||p and < ||/(2n-)||p. 

Therefore it remains to bound the contribution from / = {1, . . . , A;} (in the case 1 = 0). But in this case 
r = 0, so \Vr^{nicfj)\ = \Ujcfj\ = |(99^/J)| < C||/||^ and exp(-Apt(n - 2|/|)) < 1, which easily gives 
the desired estimate. 

□ 

In an analogous way, replacing Fact 15.91 by Fact 15.101 one proves 

Corollary 5.13. Let {Xt}^^Q be the OU branching system with Xp = 2/i. There exist constants C,c such 
that for any canonical f G 'P(M"'^) and r G (0, 1) we have for t > 1, 

E(lEMt\Xt\r''/'\Unm{lX,\>re^r'^}^^ 

We also have an analogous statement in the supercritical case. 
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Corollary 5.14. Let {Xt}^^Q be the OU branching system with Xp > 2/x. There exist constants C,c such 
that for any canonical f G P(M"'^) we have 



E (e-"("''-'^)*|[/r(/)|) < Cexp{c||x||} \\f{2n.) 
Proof. Using the notation from the proof of Fact 15.121 we get 



/C{Z+l,...,fc} 



5 E 

/C{/+l,...,fc} 



=Ap|/|t 



i(Ap-At)t+Ap|/|tTA|-f1 



Let us note that by Ap > 2/i and (j4ip we get n(Ap — /i) — Ap|I| > ri(Ap — /i) — 2(Ap — > j/'^KAp — /i). 
Thus the summands on the right hand side above for / 7^ { 1 , . . . , /c} can be bounded using CoroUary 15.71 
and Fact 15.111 by 

Cexp(c||x||)||n,c/j||^ < C2exp(c||x||)||/(2n-)|| 

(the last inequality is analogous as in the proof of Fact I5.12p . 

The contribution from / = {1, . . . (in the case / = 0) also can be bounded like in Fact 15.121 
Namely ^ = 0, so \vP{Iiicfj)\ = |n,c/j| = \{^\fj)\ < and exp(-n(Ap - + Ap|/|t) < 

exp(— (n — 2|/|)(Ap — fi)) < 1, which easily gives the desired estimate. □ 

Let Hi,fi2 be two probability measures on M, and Lip(l) be the space of 1-Lipschitz functions M 1— ?• 
[-1,1]. We define 

m{fii,ij,2) := sup \{g,fii)-{g,i22)\. (42) 
ggLip(l) 

It is well known that m is a distance metrizing the weak convergence (see e.g. [9l Theorem 1L3.3]). One 
easily checks that when f^i, fi2 correspond to two random variables Xi,X2 on the same probability space 
then we have 

mifii,fi2) < ||Xi -X2II1 < ^J\\X^-X2\\2. (43) 

Fact 5.15. Let {Xt}^yQ be the OU branching system starting from x and Ap < 2/i. For any n>2 there 
exists a function In ■ 1R+ 1— )■ M+, fulfilling lims_5.o ln{s) = and such that for any canonical /i, /2 G 'P(M"'^) 
and any t > 1 we have 

m{fii,fi2) < /n(||/i(2n-) - /2(2n-)||p), 

where ~ \Xt\~^^'^U]^{fi), fJ.2 ~ \^t\~"'^'^U^{f2) (i^^ U -statistics are considered here conditionally on 
Ext""). 

Proof. Let us fix (7 G Lip(l). We consider 
L := I (5, - {g,^i2) I = % (\Xt\-^I^U^{f^)) - ¥.g (\Xt\-^I^Unf2 



< E 



g (ixr"/'c/r(/i)) - 9 {\xtr/^unf2)) 
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Let h{x) := /i(2nx) — f2{2nx), take r := and assume that r < 1. Then by Coronary 15.121 we get 

E \g (iXtr^/^Unh)) - 9 (|Xr"/'f/r(/2)) I l{|x,|>re-^.*} < Cn exp {c„||x||} \\h\\]l\ 
On the other hand, 

\g {\Xt\-^/^Unh)) - g {\Xt\--'^U^{h)) \ l{|x,|<re-.*} < 2||5llooF (|Xi|e-V < . 
Since on Exf^ we have \Xt\ > 1 and \Xt\e~^p^ converges to an absolutely continuous random variable, 

lim supP (iXAe'^p^ < r) = 0, 

r->-0+ t V / 

which ends the proof. □ 

An analogous proof using Corollary 15.131 gives a counterpart of the above fact in the critical case. 

Fact 5.16. Let {Xt}^^^ be the OU branching system starting from x and Xp = 2^. For any n>2 there 
exists a function In '■ IR+ ^ fulfilling lim^^o ln{s) = and such that for any canonical /i, /2 G ViW^"^) 
and any t > we have 

m{fii,fi2) < U\\fi{2n-) - /2(2n-)||p), 

where fii ~ (t|Xt|)~"/^C/"(/i), /.i2 ~ |)~"/^i7"(/2) (the U -statistics are considered here conditionally 
on Ext^). 

5.3 Proofs of main theorems 

5.3.1 The law of large numbers 

Proof of Theorem \3.1[ Consider the random probability measure fit = \^t\~^Xt (recall that formally we 
identify Xf with the corresponding counting measure). By Theorem 14.11 with probability one (conditionally 
on Exf^), fit converges weakly to ip. Thus, by Theorem 3.2 in [3], /xf" converges weakly to c/j®". But 
(/, ") = \Xt\~'^Vp{f), which gives the almost sure convergence \Xt\~^V-P{f) — )■ (/, tp). Now it is enough 
to note that the number of "off-diagonal" terms in the sum (f39|) defining f/"(/) is of order and use 

the fact that \Xt\ — )• oo a.s. on Exf^. 

The proof for / E "P (M"'^) follows directly from the central limit theorems from Section [6] (we will not 
use Theorem [XT] in their proof). □ 

5.3.2 CLT — slow branching rate 

Proof of Fact The sum (jl2p is finite hence it is enough to prove our claim for one L{f,j). Without 
loss of generality let us assume that E^ = {(1, 2), (3, 4), {2k — 1, 2k)} and Ay = {2k + 1, . . . , n} (we recall 
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notation in Section [2.ip . Using the same notation as in (llip we write 



J{z2k+l, ■ ■ ■ ,Zn) 



,U2,.. .,Un) 




H{f){zi,Zi,Z2,Z2, ...,Zk, Zk, Z2k+1 Zn)tl2 (d^i) /i2 (dZ2) • • • /-f2 (dZfc) • 



D JD JD 



where D := M+ x M"'. We know that L{f, 7) = In-2k{J{x2k+i, ■ ■ ■ , Xn))- By the properties of the multiple 
stochastic integral [171 Theorem 7.26] we know that EL(/, 7)^ < ^niG{2fc+i n} Id /^il^-^i)) \ J{z2k+i, ■ ■ ■ , ^n)P- 
Therefore we need to estimate 

H{f)iz\, z\, z\, z\, . . . , z\, z\, Z2k+1 ■ ■ ■ , Zn)H{f){zf, Z^, Z2 i • • • ; Z^, Z^, Z2k+l • • • , 2;n) 

/DJD JD 

fj.2 {dzl) fi2 {dzl) ... ^2 (d4) (dzi) ^2 (d^a) . . . /i2 (d^^) /ii (dz2fc+i) ■ ■ ■ /^i (dzn) < (*). (44) 

We will now estimate H{f){zi, Z2, ■ ■ ■ , Zn)- Let Yi,Y2, . . . , y„ be i.i.d., 1^ ~ 99. We define Yi{t) = x^e"'^* + 
ou{t)Yi. Recall our notation Zi = {si,Xi) and let / = {i € {1, . . . , n} : Sj > 1}, = {1, . . . , n} \ /. Let 
/ := // be defined according to (I28p . 

Using (j23p . Lemma [5. 31 the assumption that / is canonical and the semigroup property, we can rewrite 
(EOD as 



H{f)izi,Z2, ...,Zn):= EfiYiisi), y„(s„)) = Ef{iYi{s^ - iY,isi))ieic) 

= E Yl {-l)^^^'''f((Y,{s,-l) + e,{Y,-Yi{si-l)) 
{e,)isie{o,i}^ 



i&i \ / iei^ 



Yi 



By Lemma 15.21 we have (in order to simplify the notation we calculate for d = 1, the general case is an 
easy modification) 



\Hif){zi,Z2,...,z„)\<\\fyEll 



Y, 



< 



V 



Y, 



gS/» + e~yi) dyi 



\Y^is^)\ 



< II/IIp (n (^1 ^My^{si - 1)) - exp(y,)l + El eM-Yiisi - 1)) - exp(-yi)|) V J] e 

Obviously for any x, y € M we have max (exp(x), exp(y)) < exp(x) + exp(y). Therefore by the mean value 
theorem we get 



E| exp(yi(si - 1)) - exp(yi)| < E|y,(si - 1) - y^l exp(yi(si - 1)) + E|yi(si - 1) - y^l exp(yi). 
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Using the Schwarz inequality and performing easy calculations we get 



E|exp(yi(s,- 1)) -exp(yi)| 

< VE(yi(si - 1) - Yif (VlEexp(2yi(s, - 1)) + VlEexp(2yi)) < exp (2|xi|) e-^^\ 

Similarly 

E| exp(-yi(si - 1)) - exp(-yi)l < exp {2\xi\) e-^'\ 
Since we also have Eexp(|yj(sj)|) < exp(|a;j|), we have thus proved that 

Z2, . . . , < ll/ll^e^ELi J] e-'^^' (45) 

is/ 

We use the above inequality to estimate (|44p . To this end let us denote Dq = [0, 1) x M"', Di = [1, oo) x 
M^. To simplify the notation let us introduce the following convention. For subsets /i,/2 C {1, . . . , fc}, 
/s C {2A; + 1, . . . , n} and i E {1, . . . , n} we will write = 1 if i G /j and = otherwise. Let us 
also denote 

F(z) = R{j){z\,z\,z\,z\,.. . , zl, zl, Z2k+1 ■ ■ . ,Zn)H{f){zl,zl,zl,zl, . . . , zl, zl, Z2k+1 ■■■,Zn) 
for Z= {zl, . . . , zl, zf, . . . , zl, Z2k+1, ■ ■ ■ , Zn)- 

With this notation we can estimate (14411 as follows 



(*)= E E (n/ A^2(dz.^))(n / /^2(dz,2)V f[ [ /.2(dzo)i^(z) 

=: E E ^ai>^2,/3). 

/l,72C{l,...,fc}/3C{2fc+l,...,n} 

Now, using (j45p in combination with the Fubini theorem, the definition of the measures /ij (given in 
Section [3T]) and our assumption Ap < 2^, we get 



To conclude the proof we use the fact that / i— t- L{f,-y) is linear and || • ||p is a norm. □ 

Proof of Theorem \3.SX For simplicity we concentrate on the third coordinate. The joint convergence can 
be easily obtained by a modification of the arguments below (using the joint convergence in Theorem 14.21 
for n = 1). In the whole proof we work conditionally on the set of non-extinction Exf^. 
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Let us consider bounded continuous functions /{, /g, . . . , / = l,...,m, which are centred with 
respect to 99 and set // := ^^Li/i ai^d / = "^i^i fi- In this case the [/-statistic ([2]) writes as 

m \Xt\ 

UtU) = E E f{{XS,))fi{XS2)) . . . fn{Xt{ln)). 

1 = 1 il,i2,...,in = l, 

ij^ik, for jj^k 

Let 7 be a Feynman diagram labeled by {1, 2, . . . , n}, with edges and unpaired vertices A^. Let 

Sh):= E f{Xt{ii),Xt{i2),...,Xt{in)). 

ij=ik, if {j,k)GE-, 



Decomposition (i27l) writes here as 



(46) 



where the sum spans over all Feynman diagrams labeled by {1, 2, . . . , n}, and the remainder R is the sum 
of ^-statistics corresponding to partitions of {1, . . . ,n} containing at least one set with more than two 
elements. We will prove that _^ q_ To this end let us consider (\J^. Ar) U (U,. Br) U (U^ Cr) = 

{1,2,... ,n}, a partition in which \Ar\ > 3, \Br\ = 2 and \Cr\ = 1- Each set in this partition denotes 
which indices are the same (e.g. if {1,2,3} = Ai then ii = 12 = is in our sum). By our assumption there 



is at least one set Ar- For any / = 1, . . . , m, we write 



\X^\~\AA/2 ^ fUXtiikXfLiXtiik,)) ■ ■ ■ fijXtitU) 



\ 



I 



( 



n 



n 



i^tr'E/fci(^*(^fci))4(^*fe)) 

r 

y where{fci}=Cr 

By Theorem 14. II the first term converges to and the second one converges to a finite limit. The third 
term, by Theorem 14. 2| converges (in law) to the product of Gaussian random variables. We conclude 
that K\ —7''^ 0. Thus only the first summand of (|46|) is relevant for the asymptotics of [/"(/). 
Consider now 



i^tr("/^)5(7) = E n \M~'Y.fm^)fk{Xt(i)) n i^r^/^E/'(^*») • 



1=1 {j,k)eE.y 



1=1 



1=1 



Let us denote ZriXt) := l^tT^/^ EEi /^(^tlO)- By Theorem and the Cr amer-Wold device we get 



that 



iZflit))l<j<n,l<l<m -^'^ (G fl)l<j<n,l<l<m, 
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where {G j-i)i<j<n,i<l<m is a centred Gaussian vector with the covariances 

Cov(G^,,G^,) = (<^,/j/i)+2Ap (95, (e(V2).r./j) (e(V2)«r./i) ) ds. (47) 



On the other hand by Theorem 14.11 one could easily obtain 

\Xt\-'^f,iMi))fi{Xti{})^{fjfi,ip), a.s. 

i=l 

We conclude that 

By decomposition ()46p and the considerations above we obtain 

|Xi|-("/2)*C/«(/) ^(-1)K7)^(^) 

7 

We will now show that L is equal to Li{f) given by (jl2p . By linearity of Li{f) it is enough to consider 
the case of m = 1. We will therefore drop the superscript and write /j instead of fl- We denote by 
PifiJj) ■= Jo^ifi fj)iz,z)fJ-2idz). By gZD and definition of /i2 given in ([II]) 

We now adopt the notation that r/ C 7 when E"^ C -E-y and write 

7 VC-y ij,k)eE^\En {j,k)£Er, reAj 

=E n mjk)Yi(-^y^'^"^^'^ n lEG^.G/. n%- 

{j,k)eE^ l^ri {j,k)&E^\En r&A-y 

Let us notice that the inner sum can be written as 

where cr runs over all Feynman diagrams on the set of vertices A^. Thus by [17^ Theorem 3.4 and Theorem 
7.26] this equals I\Ar,\ {H{(^i^Anfi)) and in consequence 

L = Y. n P{f,Jk)I\A,^{H{®^eAJi)). 
V {j,k)(iEr, 

It is easy to notice that in the case of / = ^f^ifi the expression above is equivalent to (fT2]) . Let 
us now consider a function f £ V. We put h{x) := f{2nx). By Lemma l5.ll we may find a sequence of 
functions {/ifc}^. C span{A) such that ^ h m. V . Next we define /^(rE) := h\^{xl2n). Now by Fact 15.161 
we may approximate \Xt\-^'^l'^^U'^{j) with \Xt\-^'^^'^^Ul\h) uniformly in t. This together with Fact E^] 
and standard metric-theoretic considerations concludes the proof. □ 
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5.3.3 CLT — critical branching rate 

Proof of Fact l374\ and Theorem \3.5[ As in the subcritical case, we will focus on the third coordinate. 

Recall the notation of Lemma l5.ll and note that A C where denotes the algebraic tensor 

product (or more precisely its standard realization subspace of L2iW"^,U7=i^idxi))). Thus the 

operator L2 is clearly well defined on ^ n Can. In the course of the proof of the limit theorem we will 
show that L2 is bounded on this space equipped with the norm || • ||-p and thus by Lemma [5.11 it extends 
uniquely to a bounded operator L2 : Can — )■ L2{^,J-,P)- 

The proof is slightly easier than the one of Theorem 13.31 as, because of larger normalization, the 
notion of [/-statistics and ^-statistics coincide in the limit. Let us consider bounded continuous functions 
/|, /2, . . . , fn, / = 1, . . . , m, which are centred with respect to if and denote / := X^JI^i ®F=i/i • By ()46|) 
we have 

unf)-vrif)= Yl {-iy^''^s{^) + R, 

7 with at least one edge 

simply by the fact that the Feynman diagram without edges corresponds to V^{f). Analogously as in 
the proof of Theorem 13.31 we have i? — )■ 0. Let us now fix some diagram 7. Without loss of generality we 
assume that Ey = {(1, 2), (2, 3), . . . , {2k - 1, 2A:)} for > 1. We have 

m 

{t\x,\r-/'s{^) = Y^]j[{\x,\t)-'{xuLfL+i)) n {{\xt\tr'/'{x,ji)). 

1=1 i<k i>2k 

By Theorem 14.31 each of the factors for i > k converges in distribution, whereas by Theorem 13.11 each 
factor for i < k converges in probability to 0, in consequence so does {t\Xt\y^'^S{j), which shows that 
(i|Xi|)-"/2c/f (/) and {t\Xt\y/'^Vl'{f) are asymptotically equivalent. 

Let us denote Zfi{t) := ^JJ' f^.{Xt{i)). By Theorem O and the Cramer- Wold device we 

get that 

^^^\Xt\,{Z ji{t))i<i<n,l<rnj '^'^ {G ji)i<i<n,l<rnj , 

where W has exponential distribution with parameter [2p — l)/p, {G fi)i<_n i<m is centered Gaussian with 
the covariances given by (jl4p . moreover W is independent of {G fi)i<j<n,i<m- Thus 

Jj _ _ 5 _ 

(m n \ 
^'En^/] j ={w,L2{f)). 

In particular (te^f*)-"/Vt"(/) VF"/2^2(/). Moreover, by FactEIOl ||(te^''*)""/Vt"(/)||2 < G\\f\\v 
for some constant G and since this bound is independent of t we get ||VF"''^L2(/)||i < C||/||-p. Recall 
that for any n there exists a constant c > such that for any > and any A^-variate polynomial 
P{Gi, . . . , Gn) of degree n in a Gaussian random vector (Gi, . . . , Gn), we have 

P(|P(Gi,...,G^)| >c||P(Gi,...,G;v)||2) >c. (48) 
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This is a well known fact in the theory of polynomial chaos and follows e.g. from a combination of Theorem 
3.2.5, Theorem 3.2.10 and Proposition 3.3.1. (i.e. the Paley-Zygmund inequality) in [?]• Note that c does 
not depend on N, which is crucial in our application. Let f £ Af] Can. Now, for sufficiently large D, by 
Markov's inequality and the fact that W is with probability one strictly positive, 

n\L2{f)\ > D^fWv) < ¥{\W^/'L2{f)\ > DWfWr) + P(W^-"/2 > d) < ^ + ¥{W~^/^ > D) < c, 

which, together with (I48p . implies that ||L2(/)||2 ^ and in consequence the operator L2 is 

bounded on A Can. In particular it admits a unique extension to a bounded operator L2 ■ Can — )■ 
L2in,T,F), which proves Fact [Ml 

Theorem 13.51 follows now by the already established case of / G span{A), Fact 15.151 and standard 
approximation arguments. 

□ 

5.3.4 CLT — supercritical branching rate 

Proof of Theorem \3.8[ Again we concentrate on the third coordinate. The joint convergence can be easily 
obtained by a modification of the arguments below (using the joint convergence in Theorem 14. 41 for n = 1). 

First, note that [/-statistics and ^-statistics are asymptotically equivalent. The argument is analogous 
to the one presented in the proof of Theorem 13. 5| since under assumption Ap > 2/i we have 

\Xt\ 

-— r-^ — )■ a.s. 

exp(2(Ap - fi)t) 

as n —7- 00 and consequently we can disregard the sum over all multi-indices (ii, . . . in which the 
coordinates are not pairwise distinct. 

Let us consider bounded continuous functions fi, f2, ■ ■ ■ , fn - IR'^ — ?• IR, I = 1, . . . , m, which are centred 
with respect to ip and denote / := X][2=i ®^=ifi- By Theorem 14.41 for n = 1 we have 

m n m n 

^-niX-,)tyn^f^ ^ ^ JJ (^g-n(A,-M)ty^i(j/)) ^ ^ J] ( (g^^^ /i , y,) o H^j = L,{f), in probability. 

1=1 i=l 1=1 1=1 

Before our final step we recall that the convergence in probability can be metrised by d{X,Y) : = 
^ (]x^T7^+t) — ^1"^ ~ Y\'^. Let us now consider a function / G "P. By Lemma 15.11 we may find a 
sequence of functions {fk} C span(A) such that — t- / in "P. Now by Corollary 15.141 we may approxi- 
mate e~"^^'P~'^^^Ul^{f) with e~"'^'^p~^'>^U]}{fk) uniformly in t in the sense of metric d. Moreover, one can 
easily show that limfc^^_oo c^(-['3(/fc), -^3(/)) = 0. This concludes the proof. □ 

6 Remarks on the non- degenerate case 

Let us remark that as in the case of J7-statistics of i.i.d. random variables, by combining the results 
for completely degenerate [/-statistics with the Hoeffding decomposition, we can obtain limit theorems 
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for general [/-statistics, with normalization, which depends on the order of degeneracy of the kernel. 
For instance, in the slow branching case Theorem 13.31 the Hoeffding decomposition and the fact that 
Ilfc: 'P(M'"^) P(M"'') is continuous, give the foUowing 

Corollary 6.1. Let Xp < 2fi and f G 'P(]R""') be symmetric and degenerate of order k — 1. Then condi- 
tionally on Ext", |Xt|-("-'^V2)[/n(j _ converges in distribution to (^)Li(nfc/). 

Similar results can be derived in the remaining two cases. Using the fact that on the set of non- 
extinction \Xt\ grows exponentially in t, we obtain 

Corollary 6.2. Let Xp = 2fi and f G P(]R"'^) be symmetric and degenerate of order k — 1. Then condi- 
tionally on Ext", t-^/'^\Xt\-^''-^/'^^Ul'{f - converges m distribution to (^)L2(nfc/). 

Corollary 6.3. Let Xp > 2fi and f € V{W^'^) be symmetric and degenerate of order k — 1. Then condi- 
tionally on Ext", exp(— (Apfi — fik)t)U]}{f — (/, 9^®")) converges in probability to (^) VF"~'^'L3(nfc/). 
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